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We develop an alternative approach to time independent perturbation theory in non-relativistic 
quantum mechanics. The method developed has the advantage to provide in one operation the 
correction to the energy and to the wave function, additionally we can analyze the time evolution 
of the system. To verify our results, we apply our method to the harmonic oscillator perturbed by 
a quadratic potential. An alternative form of the Dyson series, in matrix form instead of integral 
form, is also obtained. 
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I. INTRODUCTION 



The Schrodinger equation [l| is the main equation in non-relativistic quantum mechanics. Although it has been 
widely studied since its introduction, there are only few cases in which it can be solved exactly. Typical examples 
of potentials where an exact analytical solution is known, are the infinite well, the harmonic oscillator, the hydrogen 
atom 0, 0| and the Morse potential The great majority of the problems related with the Schrodinger equation 
are very complex and can not be solved exactly, as is the case of the cosine potential [5| for instance. When an exact 
I analytical solution can not be found, we are forced to apply approximation methods [3, 6], that when correctly used, 
i*- give us a very good understanding of the behavior of the quantum system. For the sake of clarity, we will briefly revise 
one of those methods: the time independent perturbation theory, also known as Rayleigh-Schrodinger perturbation 
theory. This method has its roots in the works of Rayleigh and Schrodinger, but the mathematical foundations were 
only set by Rellich in the late thirties of the past century (see Simon p} and references there in). This method has 
Q" 1 been applied with great success to solve a vast variety of problems such that, through its continuous implementation, 
a lot of techniques have been developed, which go from numerical methods [gj to those more mathematical and 
I | fundamental, as convergency problems [tl fioj . The Rayleigh-Schrodinger theory is appropriate when we have a time 
. independent Hamiltonian, that can be separated in two parts, as follows: 

(N ! H = H Q + H pi (1) 

r- ■ 

where Hq is the so-called non-perturbed Hamiltonian, and it is usually assumed to have known solutions, this is, its 
c"^i • eigenvalues 

ff |n<°>) = 4°V (0) ) (2) 



are known. The second part of the Hamiltonian, H p , is small compared to Hq; thus H p is called the "perturbation", 
because its effect in the energy spectrum and in the eigenfunctions will be small. To be more explicit, it is usual to 
write H p in terms of a dimensionless real parameter A, which is considered very small compared to one 

H p = XV (A<1). (3) 

Here, we propose an alternative perturbation method based on the evolution operator for Hamiltonian (JlJ, i.e., the 
exponential of the complete Hamiltonian. Making use of the Taylor series of the exponential and operator techniques 
we will cast the solution into a form that has powers of tridiagonal matrices in it. By using as an example a perturbative 
potential for which we know the solutions, we will be able to compare our method with the Rayleigh-Schrodinger 
perturbation theory. Our method will allow also to cast the Dyson series, which are written in multiple integral forms, 
as a series of powers of tridiagonal matrices. 

We will proceed then as follows, in section II, we present a brief summary of the standard time independent 
perturbation theory, emphasizing the expressions obtained for the first and second order corrections. As is well 
known, two expressions are obtained at each order, one for the energies and one for the wave functions. In Section III, 
we introduce our matrix approach to perturbation theory that will allow us to obtain a single solution that contains 
both, the energy and wave function corrections. We will introduce there all the corrections to the wave function 
that will allow us to generate later the Dyson series, and in this way give also a new expression of it but now in 
terms of a matrix series. In Section IV we compare the our method with the standard perturbation theory by solving 
the harmonic oscillator with a quadratic perturbation. Although this seems redundant, because this case has an 
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exact solution, it will allow us to compare both methods by doing an expansion of the exact solution in terms of 
the perturbation parameter, A. In Section V we apply our formalism to rewrite the Dyson series in matrix form and 
Section VI is left for conclusions. 

II. PERTURBATION THEORY 

As we already mention, we have to solve the eigenvalue problem given by ([T]). Standard perturbation theory 
produces the following expressions for the eigenvalues 

£;„ = 4°)+AA( 1 ) + A 2 A( 2 ) + ..., (4) 

and the eigenfunctions 

|n) = |n (0) ) + A|n«) + A 2 |n (2) ) + • • • . (5) 

In the two previous equations the super index (j) indicates the correction order. The first and second order corrections 
for the energy are 

A« = (n^\V\n^)=V nn (6) 

k^n &n - &k 

where we have defined 

V kn = {kW\V\nM). (8) 
For the wave functions the first two order corrections are given by the expressions 



k^n ^k 

and 



xv _ V \k^)V kn 
I 2^ ^(0) _ E (o) W 



(2)\ -TV \k (0) )VkmV mn _ ^ |fcW)KnVfen nm 

2^ 2^ ,„(o) p (o)w p (o) p (oh 2^ rp (o) p (0)s 2 1 ' 

It is clear that the previous expressions can be used in the non-degenerate case, where we always have En ^ ^ 
In the degenerate different treatment is needed. 

III. MATRIX APPROACH TO THE PERTURBATION THEORY 
A. First order correction 

We now find an approximate solution for the complete (time dependent) Schrodinger equation; i.e., we will solve 
approximately the equation 

i^)=H\^)=(H + \v)\i>). (11) 

The formal solution may be written as 

W)) = e - lt(do+M>) |V(0)>, (12) 

that by expanding the exponential in a Taylor series 

n=0 
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is obtained. If we develop the binomials inside the summation, rearrange terms and cut the series to first order in A, 
we have 

w*)> - (E^^ + ^E^E^r 1 '^) hmo)>- (14) 

\n=0 ' 71=1 ' fc=0 / 

The key ingredient of the method we introduce in this contribution is the matrix 

It is very easy to get convinced that the following relations are satisfied 

(M)i, 2 = V, 
{M 2 ) ia = H V + VH , 



n-l 



k=0 



Thus, equation (TTJ| can be written as 



W)) 



E^w+aE 



(-it)»(M») (li2) 



,n=Q 



n=l 



I0(o)). 



(16) 



As M° = J (where / is the [2 x 2] identity matrix), we trivially have that 

(M°) (1)2) = 0, 
and substituting this term in equation (1161) we obtain 



\m) « + A ( e - a/t ) (1 , 2) ) 10(0)) = e-^^(0)) + A (e- lM *) (1 . 2) I0W) 



(1,2) 



(17) 



Note that in the last expression we have separated the approximate solution in two parts; the first part is the solution 
of the non-perturbed system, that is well known, and the second part is the first order correction to the wave function. 
We now show how the correction to first order may be calculated, for this we rewrite equation (1171) as 



|V(*)>»l^ 0) (*)>+A(||^»i, 
where we have defined the matrix wave junction 



II0 P ) 



( I -0(1,1)) 10(1,2)) 
V 10(2,1)) 1^(2,2)) 



Deriving (|17[) and (|18l) with respect to time and equating them, we arrive to the equation 



d 



d 



i^l0 (O) (*)> + *A-(||V p ))i,2 = H e- iHat \m) + A (Me- UM I\i>(0))) h 



Equating powers of A, we have 



ig- t m\ 2 = (Me- itM i\m)) 1 



We have to solve now equation ([2"T]) . or equivalently the system 



; 9_ ( 10(1,1)) 10(1,2)) ^ = Me -UM ( 10(0)) 



8t \ 10(2,1)) 10(2,2)) 



|0(O)> 



(18) 



(19) 



(20) 



(21) 



(22) 
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Integrating this equation, we have 



to finally write the differential equation 



with the initial condition 



m _ p -um ( 1^(0)) 



i-||^>=M||^}, 



The needed solution is associated with the second column of matrix ||?/> p ), thus we write 



(23) 



(24) 



(25) 



<L ( 1^(1,2)} 



M 



1-0(1,2)) 
I -0(2,2)) 



at V 1^(2,2)) 

but, because M is a tridiagonal matrix, ithe system may be directly integrated. We show it by writing explicitly 

d 



«^|-0(1,2)) = H \lp(l,2)) + V|^( 2 ,2)>, 



and 



with the initial condition 



^1^(2,2)) = i^o|-0(2,2)), 



( 

V I0(o)> 



Because we know the solution for H , equation (j2"7) is solved trivially, 



I V>(2, 2 )> =e- im °\mh 



(26) 



(27) 



(28) 



(29) 



that substituted in (|26[) . allows us to write 



d 



^IV>(i, 2) ) = ffo|V(i,a)> + Ve-^im) 



Making the transformation |V>(i,2)) = e °\ < Pi,2{ x ))i we arrive to the equation 



that can be easily integrated to give 



I0i,2) = -i 



e iHotl Ve- iHotl dti 



IV»(0)>, 



and by transforming back the solution to the first correction may be obtained 



-iHnt 



e iHotl Ve- iHotl dti 



I0(O)>- 



(30) 



Up to here we have produced a first order correction for the wave function with no assumptions on Hamiltonian 
degeneracy, therefore making this first order correction valid also for degenerate Hamiltonians. 
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As the eigenfunctions of the non-perturbed Hamiltonian constitutes a complete orthonormal set, we can write the 
closure condition as / = J2k \k^°'){k^\, where |fc ( - ' ) ) are the eigenstates of the unperturbed Hamiltonian and insert 
the identity operator written in this way, in equation (|30|) to arrive to 



|^, 2 >=-^|fc (0) }<A ; ( )|e- < ^ i 



|V(0)). 



(31) 



For simplicity we set the initial condition |^(0)) = \n^). By defining 



(1) >* = E 



\k {0) )V kn ; E (o h 



E (o) _ E (o) 



and using expressions © and ([9]), we finally get 



1^1,2 



(32) 



(33) 



We note that the first order corrections to the energy and to the wave function are already contained in the last 
expression. Additionally we can show that we can write 



P ) = e -itH fl -ifodtiV(ti)\ ]m)) 



with 



iflnti f>„— iffntl 



V(h) = e lHotl Ve 



(34) 



(35) 



B. Second order correction 

We will find now a similar expression for the second order correction. Expanding again (|13[) in Taylor series and 
keeping terms to A 2 , we have 



Vn=0 



n — 1 n—k 



k=0 



n=2 k=l j=0 

In analogy with the first order correction, we define now the matrix 

H V 
\l = | H Q V 
H 



It is very easy to see that in this case 



(m 2 ; 



1,3 



V 



(M 3 ) li3 - H Q V 2 + VH V + V 2 H Q , 



(36) 



(37) 



11,3 



i—l n — k 



Y^Y^H^-^VHlViHof- 1 -, 



k=l j=0 



so we can write equation (|3o|) as 



\m) * ( e- i6 °* + A (e~ iMt ) a2) + A 2 £ (M™) 13 ] |^(0)) 



'(1,2) ' " n \ 

n=2 



(38) 
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and knowing that 



we obtain 



Inserting now the matrix 



(M°) 13 = (M) 13 = 0, 
\m) * ( e -^' + A( e -^) (12) +A 2 (e-^) (13) )|^(0)). 



10(1,1)) 1-0(1,2)) 1^(1,3)) 
! J P ) = \ 10(2,1)) 10(2,2)) 10(2,3)) 
10(3,1)) 10(3,2)) 10(3,3)) 

where the first order corrections 10(1,2)) and the second order corrections |0(i,3)) are included. 
Expanding \ip) to second order in A, we get 

I0)«|0 (O) W)+A(||0 P ))1,2 + A 2 (||0P)) 1 , 3 . 



(39) 



(40) 



(41) 



Following the same procedure as in the first order correction case, we derive equations (f31)|) and (14"T|) with respect to 
time, and equate the corresponding equations to obtain 



^I0 (O) (*))+^(||0 P )) 1)2 + A 2 |(||0^) 1)3 = H o e-^ t |0(O))+A(Me-" M 7|0(O))) 1 . 



at 



at 



A 2 (Me- itM /|0(O))) 1 



(42) 



Equating powers of A 2 , we can establish that 



d 



i^(||0 p ))i,3 = (Me- itM /|0(O))) 1 



or equivalcntly 



with the initial condition 



II0 P (O)) 



i-m = M\\r), 



10(0)) 
|0(O)) 
|0(O)) 



(43) 



(44) 



Equation (|4"4")l is similar to (pM)) , thus we can again proceed as in the first order case, choosing the third column in 
both sides of the equation and getting the differential equations system 



10(1,3)) 
10(2 
10(3,3)/ 



10(1,3) 
10(2 
10(3,3)/ 



with the initial condition 






|0(O)) 



The correction we were looking for is then 

101,3) 



-ie 





r - z 1 * 1 


[ e iAotl V 


_ ie -iH 0tl / 


Jo 


Jo 



e iHot2 Ve- i " ot2 dt 2 \ij}(0)} 



dt 1 . 



(45) 



(46) 



(47) 



In this equation we note that the expression inside the square brackets is the first order correction, then using (|3"3")l 
we can write 



101,3) 



iH t / JHoti 



V 



o 



^ (|nW)-^|^))AW-|nW) tl ) 



dti. 



(48) 
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Inserting the identity operator / as we had defined it in the first order correction, we have 



with 



and 



k 

,-i£< 0) t 



IVi.s) = -*V|/fc (0) )(fc (0) |e- 4Hot / e iHotl V 
^ Jo 

\k^)V km V r 



(| n (i))_ iil | n (o) )A (i))_| n (i) )ti 



dt-y, 



E E ,^01 ,0)10, ,0i- ^ F ,u, M 



z* 2 (A«) 2 - ]T 



it|A;W)Ak 1 V fcn 
p(o) p (0) 



|n< a >)t + |<% 



(2)\ 



EE 



|fc (0) )Vfc m Kn 



E 



p(0) _ p(0)x 2 



e 



(49) 



(50) 



EE 



|jfc(°)> 



e fe Vfcm Vm 



F (0) _ F (o) 



e «(4 0) -^ 0) ) dtl 



(51) 



In the same way as in the first order correction case, equation (|40[) gives us the correction to the wave function and 
to the energy in one operation; however, it generates the additional expressions \n^ 2 \t)) y |n^(t)). 
We can show that it may be written in the compact form 



1 -if*dhV(ti) H) 2 Jq^i/o 1 dt2V(h)V(t2) 



-iJ'dhVih) 



1-0(0)). 



(52) 



C. Higher order corrections 



In the previous sections, we have found the corrections to first and second order, and we have shown that the results 
are not only consistent with the traditional perturbation theory, but new terms arise. These new terms will be shown 
correct in Section IV where we compare both solutions. In rest of this section we generalize our method to higher 
orders. We propose the perturbation series 



-itM\ 



10(0)) 



(53) 



with 



ijP) 





(54) 



where m is the correction order, and 



M 



f H V 
H 

V ••• 



\ 





Ho) 



(55) 



Following the method for the first and second order corrections, we deduce the following system of differential equations 



or 



»^(||^»i,™+i = M(||^)) lim+1 , 



(56) 




(57) 



with the initial condition 




so the solution we are looking for is 



(58) 





r r u 








|^i,n+i) = -ie- lHat f e lHatl V 


-ie~ lHotl / e lHot2 V 


_ ie -iH„t 2 / . 


■dt 3 


dt 2 \ip(0)) 


Jo 


Jo 


Jo 







dt x - 



(59) 



IV. COMPARISON OF THE STANDARD AND THE MATRIX METHODS 

We will treat now the case of a harmonic oscillator perturbed by a quadratic potential. More than as an example, 
what we pretend in this section is to show that in this case, where an exact analytic solution is known, our method 
gives a better solution than the standard Rayleigh-Schrodinger perturbation theory. The non-perturbed Hamiltonian 
is 



Ho = \{p 2 +oj 2 x 2 ), 



and the perturbation potential is 



= T 



for simplicity we consider a unity mass oscillator 
The total Hamiltonian is 



(60) 



(61) 



(62) 



with 



uj = cjvT+A. 



(63) 



Physically this Hamiltonian can represent a one mode cavity, to which the oscillation frequency can be changed by 
changing the perturbation parameter. The hnal result of a modification in the initial frequency produces squeezed 
coherent states as it is described by Dutra 

First we obtain the first order correction with the matrix method we have just introduced. We do so by substituting 
the perturbation potential in equation p0[) and obtain 



l^i 



-iH t 



iH ti £2 e -iH ti 



IV>(0)>. 



By defining the position operator in terms of annihilation and creation operators, a and 



a + a 



(64) 



(65) 



we can write equation (|64|) as 



1^1,2) = -ije 



-iHnt 



IV'(o)). 



Using the relations 
and 

equation (|66p is transformed in 
^1,2) = -i-^e 



e iH t &e -iH t = fe-iut 



e iH ttf e -iH t = a t e *«* 



-iHnt 



-iHnt 



IV>(0)>, 



( e - 2 -* _ 1) + M- ( e 2 -* _ 1) + t (aat + a t a) 



t^2 



-HbJ 



2iuj 



(66) 

(67) 
(68) 



|V(0)). 



(69) 



By writing the unperturbed Hamiltonian in terms of the number operator n = a'a and using |"0(O)) 
obtain 



,(0) 



n^> ) = we 



iEi°h 



V^EH |„(o) _ 2(0)) _ >/(" + !)(" + 2) |„(o) + 2 (0) } ) _ | n (0) } ^ (2n + 



l)t 



5(0) ty /n{n- 1) 



(70) 



Comparing with (|33[) , we have that 



w = >(»-!) (0) _ 2(0)) _ V(n + !)(« + 2) (0) + 2(0) 



(71) 



A« = ^(2n + 1), 



and 



| n (l)^ t = c -iE?Lt y/n(n- 1) (0) _ ^ _ ( ,i,"".,> V .. n»-r - ■ , 



(72) 



(73) 



with #i 0) = io(n + 1/2). 

As we can see (|7Tj) and ff72]) are the corresponding expressions to the corrections to the wave function and to the 
energy to first order that are obtained when the stationary Schrodinger equation is considered, while (|73p is similar 
to (|71j) . but with the difference that this expression has time dependent factors. 



A. Exact solution 



We may see already differences between both methods. In order to compare them, we need to use an example that 
may be solved, such that, once we have the exact solution, this can be expanded on a (perturbation) parameter and 
see if it matches the results we obtained. 

In this case, the exact formal solution is given by 



m)) = e- iH * t \m)- 



(74) 



where is the complete Hamiltonian defined in (|62 
As usual, we introduce the ladder operators 



it 



—a; — i — == 
2 



and 



in terms of which the perturbed Hamiltonian (|62|) can be written as 



//,. = - f iU + ~ 



it 



It is easy to show that 



and 



Defining the squeeze operator [Tl| 



it can be shown that 



and 



Thus, the formal solution (|74p can be written as 



uj a + fit /u7(fi — fit) 



w a + fit /u7(fi — fit) 
2 JTT + \2 



i = .Sfi^t 



it = Stf&. 



W)) = e 



_ -itu>(A t A+l/2) 



MO)), 



= e- ita,S ^" +1 W 3i \il>(0)), 
= 5e-^( at<i+1 / 2 )5t^(0)), 

and multiplying by e^M^ «+i/2) e «a(a f 6+1/2) ; we get 

^(t)) = e ' ,(li ( sti + 1 /2) e *(n , n+l/2)5' e -"*'i(« t n+l/2) j 5tU(Q)| i 

We define now the time dependent squeezing operator 



S~ (A t) e 'itu(a^ d+1/2) itcD(a' a+1/2) 

e it < Ii(a t a+l/2) e i ln(l+A)(a 2 -(a t ) 2 ) e -it < i(a t a+l/2) 

i ln(l + A) ( e it "( ata + 1 /2) 6 2 e - 4 ti(a t a+l/2) 



exp 



itLD(a+a+l/2)^t^2 e -ittr)(a t a+l/2)^ 
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and we have 



taking ([74| the final form 



Su(X,t) = cxp 



ln(1 + A) [e- 2lt "a 2 -e 2 ^(a)f] 



(84) 



(85) 



Up to now we have an exact result. In order to compare with the approximation found with our method, we expand 
the operators in the above expression in Taylor series in terms of the "small" parameter A. First, we have to first 
order in A, 



Second, we write 



we keep first order terms, so 



and because 



and 



so to first order in A 



Third, we write 



and it can be shown, that 



and that 



so to first order in A 



s+(°)(o) = i, 



£W(0) 



dSHX) 
dX 



((at) 2 -a 2 ) 



A=0 



S\\) « 1 + A 
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S*(X,t) = f]^Mx n 



4 0) (o) = i 



4^(0) 



dX 



A=0 



4(A, t) « 1 + Ai[a 2 e- 2lt " - (a t ) 2 e 2lt "]. 



Finally, we arrive to the expression to order one in A 



W)) 



-itH, 



A 



(86) 



(87) 



itu) 



( 1 + g ((6 f ) 2 (l - e 2 ^) + a 2 (e- 2 -* - 1)) - A— ( h + - ) ] |^(0)> 



(89) 
(90) 

(91) 

(92) 

(93) 
(94) 

(95) 

(96) 



that when In' )) is taken as initial condition, formula (l69l) is recovered, that is the one obtained from the method 
introduced here. 

The development to second order is long and bothersome and we do not present it here, but we recover the expressions 
found with our method in the previous section. 
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V. THE DYSON SERIES IN THE MATRIX METHOD 

It is well known that in terms of the Dyson series [13j , the wavefunctions of the perturbed problem are written as 



\iP(t)} = e - lMot T<exp 



Wo)>. 



(97) 



-iX / dt x V{ti) 
Jo 

where f is the time order operator [ljj]; i.e, if we have the time dependent operators A(t) and B{t) then 

f[A( tl )B( t2 )} = { f f^^if 08) 

L A{ti)B(t 2 ) if t\ > t 2 . 

On the other hand, from equation (|53|) we can write 

m)) = (e- iAot + E A " ( e ~ mt ) (i,„ + d) W°)>; (") 

so comparing equation (|97p with equation ([99ft , we derive the formula 

oo 

= e^^A«(e-*) Ml) . (100) 



T jexp 



-iX / diiV(tx) 
Jo 



which offers a matrix expansion for the Dyson operator and that links our matrix method with the Dyson series. 



VI. CONCLUSIONS 

In this work, we have developed a new technique to find approximate solutions to the Schrddinger equation. We 
used the formal solution of the time dependent Schrddinger equation (|12j) . The key ingredient is the introduction of 
the matrix M defined in (|55[) that allows us the transformation of the Taylor series for the wave function, in terms of 
products of the operators Ho and V, in a series of powers of the matrix M, that is easier to handle. 

The method allowed us to express the terms of the perturbation series in the form of integrals that depend on time 
that are restricted to the interval [0, t], as appears in (|30|) and (J4TJ) . An interesting property of these equations is 
that in the form that are presented they do not distinguish if the Hamiltonian Hq is degenerate or not, for what the 
equations that we provide for the corrections are very general expressions. 

To test our method versus the standard Rayleigh-Schrodinger theory, we solved the corrections to the first and 
second order of quadratic potential as perturbation on the harmonic oscillator, and we found that the solution not 
only contains the results of the stationary case, but we find other time dependent terms that help us to study the 
temporal evolution of the corrections. It is important to remark that the terms as (l3"2"j) , (I5U1) . (I5T1) among others could 
not be found in the conventional treatment of the perturbation theory. 

Finally, we have given an alternative expression for the Dyson series in a matrix form. 
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